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I The infrared structure of quantum gravity is explored by solving a lattice version of the Wheeler- 

o\. 

•/^ ■ DeWitt equations. In the present paper only the case of 2 + 1 dimensions is considered. The nature 

I of the wave function solutions is such that a finite correlation length emerges and naturally cuts 



off any infrared divergences. Properties of the lattice vacuum are consistent with the existence 
of an ultraviolet fixed point in G located at the origin, thus precluding the existence of a weak 
coupling perturbative phase. The correlation length exponent is determined exactly and found 
^ I to be = 6/11. The results obtained so far lend support to the claim that the Lorentzian and 

Euclidean formulations belong to the same field-theoretic universality class. 
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1 Introduction 



It is possible that the well-known ultraviolet divergences affecting the perturbative treatment of 
quantum gravity in four dimensions point to a fundamental vacuum instability of the full theory. 
If this is the case, then the correct identification of the true ground state for gravitation necessarily 
requires the introduction of a consistent nonperturbative cutoff. To this day the only known way 
to do this reliably in quantum field theory is via the lattice formulation. Nevertheless, previous 
work on lattice quantum gravity has dealt almost exclusively with the Euclidean formulation in 
d dimensions, treated via the manifestly covariant Feynman path integral method. Indeed the 
latter is very well suited for numerical integration, and many analytical and numerical results have 
been obtained over the years within this framework. However the issue of their relationship to 
the Lorentzian theory has remained largely open, at least from the point of view of a rigorous 
treatment. The main supporting arguments for the Euclidean approach come from the fact that 
the above equivalence holds true for other field theories (no exceptions are known), and from the 
fact that in gravity itself it is rigorously true to all orders in the weak field expansion. 

In this paper we will focus on the Hamiltonian approach to gravity, which assumes from the 
beginning a metric with Lorentzian signature. In order to obtain useful insights regarding the non- 
perturbative ground state, a Hamiltonian lattice formulation was introduced based on the Wheeler- 
DeWitt equation, where the quantum gravity Hamiltonian is written down in the position-space 
representation. In a previous paper [IJ a general discrete Wheeler-DeWitt equation was given 
for pure gravity, based on the simplicial lattice formulation originally developed by Regge and 
Wheeler. On the lattice the infinite-dimensional manifold of continuum geometries is replaced by 
a finite manifold of piecewise linear spaces, with solutions to the lattice equations then providing 
a suitable approximation to the continuum gravitational wave functional. The lattice equations 
were found to be explicit enough to allow the development of potentially useful practical solutions. 
As a result, a number of sample quantum gravity calculations were carried out in 2 + 1 and 3 + 1 
dimensions. These were based mainly on the strong coupling expansion and on the Rayleigh-Ritz 
variational method, the latter implemented using a set of correlated product (Slater-Jastrow) wave 
functions. 

Here, we extend the work initiated in [T] and show how exact solutions to the lattice Wheeler- 
DeWitt equations can be obtained in 2 + 1 dimensions for arbitrary values of Newton's constant 
G. The procedure we follow is to solve the lattice equations exactly for several finite regular 
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triangulations of the sphere and then extend the result to an arbitrarily large number of triangles. 
One finds that for large enough areas the exact lattice wave functional depends on geometric 
quantities only, such as the total area and the total integrated curvature (which in 2 + 1 dimensions 
is just proportional to the Euler characteristic). The regularity condition on the solutions of the 
wave equation at small areas is shown to play an essential role in constraining the form of the wave 
functional, which we eventually find to be expressible in closed form as a confluent hyper geometric 
function of the first kind. Later it will be shown that the resulting wave function allows an exact 
evaluation of a number of useful (and manifestly diffeomorphism-invariant) averages, such as the 
average area of the manifold and its fluctuation. 

From these results a number of suggestive physical results can be obtained, the first one of which 
is that the correlation length in units of the lattice spacing is found to be finite for all G > 0, and 
diverges at G = 0. Such a result can be viewed as consistent with the existence of an ultraviolet 
fixed point (or a phase transition in statistical field theory language) in G located at the origin, thus 
entirely precluding the existence of a weak coupling phase for gravity in 2 + 1 dimensions. Simple 
renormalization group arguments would then suggest that gravitational screening is not physically 
possible in 2 + 1 dimensions, and that gravitational antiscreening is the only physically realized 
option in this model. A second result that follows from our analysis is an exact determination of the 
critical correlation length exponent for gravity in 2 + 1 dimensions, which is found to be v = 6/11. 
It is known that the latter determines, through standard renormalization group arguments, the 
scale dependence of the gravitational coupling in the vicinity of the ultraviolet fixed point. 

A short outline of the paper is as follows. In Sec. 2, as a general background to the rest of 
the paper, we briefly describe the formalism of classical canonical gravity, as originally formulated 
by Arnowitt, Deser and Misner. The continuum Wheeler-DeWitt equation and its invariance 
properties are introduced as well at this stage. In Sec. 3 we introduce the lattice Wheeler-DeWitt 
equation derived in a previous paper p!], and later Sec. 4 makes more explicit various quantities 
appearing in it. This last section also discusses briefly the role of continuous lattice diffeomorphism 
invariance in the Regge framework as it applies to the present case of 2 + 1-dimensional gravity. Sec. 
5 focuses on the scaling properties of the lattice equations and various sensible choices for the lattice 
coupling constants, with the aim of giving eventually a more transparent form to the wave function 
results. Sec. 6 gives a detailed outline of the general method of solution for the lattice equations and 
then gives the explicit solution for a number of regular triangulations of the sphere. Later, a general 
form of the wave function is given that covers all the previous discrete cases and allows a subsequent 
study of the infinite volume limit. Sec. 7 focuses on one of the simplest diffeomorphism-invariant 
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averages that can be computed from the wave function, namely the average total area. A brief 
discussion follows on how the latter quantity relates to the corresponding averages computed in the 
Euclidean theory. Sec. 8 extends the calculation to the area fluctuation and shows how the critical 
exponents (anomalous dimensions) of the 2 + 1-gravity theory can be obtained from the exact wave 
function solution, using some rather straightforward scaling arguments. Sec. 9 discusses some 
simple physical implications that can be inferred from the values of the exact exponents and the 
fact that quantum gravity in 2 + 1 dimensions does not seemingly possess, in either the Euclidean or 
Lorentzian formulation, a weak coupling phase. Sec. 10 contains a summary of the results obtained 
so far. 

2 Continuum Wheeler-DeWitt Equation 

Since this paper involves the canonical quantization of gravity we begin here with a very brief 
summary of the classical canonical formalism [2] as derived by Arnowitt, Deser and Misner [3]. 
While many of the results presented in this section are rather well known, it will be useful, in view 
of later applications, to recall the main results and formulas and provide suitable references for 
expressions used later in the paper. 

The first step in developing a canonical formulation for gravity is to introduce a time-slicing 
of space-time, by introducing a sequence of spacelike hypersurfaces labeled by a continuous time 
coordinate t. The invariant distance is then written as 



ds^ = -dr^ = gf.^^dx'^dx'' = gijdx'dx^ + 2gijN'dx^dt - (N^ - g,jN'N^)dt^ , (1) 



where {i = 1, 2, 3) are coordinates on a three-dimensional manifold and r is the proper time, in 
units with c = 1. 

Indices are raised and lowered with gij{x) = 1, 2, 3), which denotes the three- metric on the 
given spacelike hypersurface, and A^(x) and A^*(x) are the lapse and shift functions, respectively. 
It is customary to mark four-dimensional quantities by the prefix ^ so that all unmarked quantities 
will refer to three dimensions (and are occasionally marked explicitly by a ^). In terms of the 
original four-dimensional metric ^g^u one has 





which then gives for the spatial metric and the lapse and shift functions 





(3) 
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For the volume element one has 

= N^, (4) 

where the latter involves the determinant of the three- metric, g = det gij. As usual g^^ denotes the 
inverse of the matrix gij . 

A transition from the classical to the quantum description of gravity is obtained by promoting 
the metric gij, the conjugate momenta tt'-', the Hamiltonian density H and the momentum density 
Hi to quantum operators, with gij and tt^^ satisfying canonical commutation relations. In particular, 
the classical constraints now select a physical vacuum state \^), such that in the source-free case 

= Hil"^) = (5) 

and in the presence of sources more generally 

f\^) = t;|^') = , (6) 

where T and Tj now include matter contributions that should be added to H and Hi. The momen- 
tum constraint involving Hi or more generally Tj, ensures that the state functional does not change 
under a transformation of coordinates x*, so that ^ depends only on the intrinsic geometry of the 
3-space. The Hamiltonian constraint is then the only remaining condition that the state functional 
must satisfy. 

As in ordinary nonrelativistic quantum mechanics, one can choose different representations for 
the canonically conjugate operators gij and vr*-^. In the functional position representation one sets 

g^j{-K) ^ 5ij(x) 7r'-'(x) ^ -ih-lQ^G- (7) 

ogij (,xj 

In this picture quantum states become wave functional of the three- metric giji'x), 

1^) ^ . (8) 

The two quantum-constraint equations in Eq. Q then become the Wheeler-DeWitt equation [H [5l 

and the momentum constraint listed below. Here Gij^^i is the inverse of the DeWitt supermetric, 
given by 

Gij,ki = ^g'^^'^ (gikgji + gugjk + agijgki) , (lo) 



5 



with parameter a = —1. The three-dimensional version of the DeWitt supermetric itself, G^^'^^x) 
is given by 

= i^(5V'+5V' + «5V) , (11) 
with parameter a in Eq. (jlOp related to a in Eq. (jlip by a = — 2a/(2 + 3a), so that a = —1 gives 
a = —2 (note that this is dimension dependent). In the position representation the diffeomorphism 
(or momentum) constraint reads 

1 2 i <7.j Vfc ^ + } ^ [9ij (x)] = , (12) 

where and Hf are possible matter contributions. In the following, we shall set both of these to 
zero, as we will focus here almost exclusively on the pure gravitational case. 

A number of basic issues need to be addressed before one can gain a full and consistent un- 
derstanding of the dynamical content of the theory (see, for example, [TJ [51 [9l [ini U] as a small 
set of representative references). These include possible problems of operator ordering, and the 
specification of a suitable Hilbert space, which entails at some point a choice for the inner product 
of wave functionals, for example in the Schrodinger form 

(*|$) = I dfi[g] ^*[g,,] <^[gij] (13) 

where dfj,[g] is some appropriate measure over the three- metric g. Note also that the continuum 
Wheeler-DeWitt equation contains, in the kinetic term, products of functional differential oper- 
ators which are evaluated at the same spatial point x. One would expect that such terms could 
produce 6^^^ (0) -type singularities when acting on the wave functional, which would then have to be 
regularized in some way. The lattice cutoff discussed below is one way to provide such an explicit 
ultraviolet regularization. 

A peculiar property of the Wheeler-DeWitt equation, which distinguishes it from the usual 
Schrodinger equation = ihdt^, is the absence of an explicit time coordinate. As a result, the 
rhs term of the Schrodinger equation is here entirely absent. The reason is of course diffeomorphism 
invariance of the underlying theory, which expresses now the fundamental quantum equations in 
terms of fields gij and not coordinates. 

3 Lattice Hamiltonian for Quantum Gravity 

In constructing a discrete Hamiltonian for gravity, one has to decide first what degrees of freedom 
one should retain on the lattice. One possibility, which is the one we choose to pursue here, is to 
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use the more economical (and geometric) Regge- Wheeler lattice discretization for gravity [12\ I13j. 
with edge lengths suitably defined on a random lattice as the primary dynamical variables. Even 
in this specific case several avenues for discretization are possible. One could discretize the theory 
from the very beginning, while it is still formulated in terms of an action, and introduce for it a 
lapse and a shift function, extrinsic and intrinsic discrete curvatures etc. Alternatively one could 
try to discretize the continuum Wheeler-DeWitt equation directly, a procedure that makes sense in 
the lattice formulation, as these equations are still given in terms of geometric objects, for which 
the Regge theory is very well suited. It is the latter approach which we will proceed to outline here. 

The starting point for the following discussion is therefore the Wheeler-DeWitt equation for 
pure gravity in the absence of matter, Eq. ([9]), 

{ - o-"!"' ^,,(x?L(^) " ^ ( " * } = ° 

and the diffeomorphism constraint of Eq. (|12p , 

I 2i5.,(x) Vfc(x) I *bi,(x)] = . (15) 

Note that these equations express a constraint on the state \^) at every x, each of the form 
i7(x) 1^) = and Hi (x)|*) = 0. 

On a simplicial lattice [Ml [I5l [TBI [13 [E] (see for example [H], and references therein, for a 
more complete discussion of the lattice formulation for gravity) one knows that deformations of the 
squared edge lengths are linearly related to deformations of the induced metric. In a given simplex 
o", take coordinates based at a vertex 0, with axes along the edges from 0. The other vertices are 
each at unit coordinate distance from (see Figs. 1, 2 and 3 as an example of this labeling for a 
triangle). In terms of these coordinates, the metric within the simplex is given by 

9vi^) = h {ll + llj-lD ■ (16) 

Note that in the following discussion only edges and volumes along the spatial direction are involved. 
It follows that one can introduce in a natural way a lattice analog of the DeWitt supermetric of 
Eq. ()lip by adhering to the following procedure [201 [21] . First one writes for the supermetric in 
edge length space 

f = G'\f) 6lf 61] , (17) 

with the quantity G^^{P) suitably defined on the space of squared edge lengths. By a straightforward 
exercise of varying the squared volume of a given simplex a in d dimensions 

V'i^) = {if det g,,{l'{a)) (18) 
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to quadratic order in the metric (on the rhs), or in the squared edge lengths belonging to that 
simplex (on the Ihs), one is led to the identification 



It should be noted that in spite of the appearance of a sum over simplices a, G^^{P) is local, since 
the sum over a only extends over those simplices which contain either the i or the j edge. 

At this point one is finally ready to write a lattice analog of the Wheeler-DeWitt equation for 
pure gravity, which reads 



with Gij{P) the inverse of the matrix G''^ (P) given above. The range of the summation over i and 
j and the appropriate expression for the scalar curvature, in this equation, are discussed below and 
made explicit in Eq. ([2T]) . 

Equations ([9]) or (f20|) express a constraint equation at each "point" in space. Here we will 
elaborate a bit more on this point. On the lattice, points in space are replaced by a set of edge 
labels i, with a few edges clustered around each vertex in a way that depends on the dimensionality 
and the local lattice coordination number. To be more specific, the first term in Eq. (|20p contains 
derivatives with respect to edges i and j connected by a matrix element Gij which is nonzero only 
if i and j are close to each other, essentially nearest neighbor. One would therefore expect that the 
first term could be represented by just a sum of edge contributions, all from within one (d — 1)- 
simplex a (a tetrahedron in three dimensions). The second term containing ^RiJ?') in Eq. ()20p is also 
local in the edge lengths: it only involves a handful of edge lengths, which enter into the definition 
of areas, volumes and angles around the point x, and follows from the fact that the local curvature 
at the original point x is completely determined by the values of the edge lengths clustered around 
i and j. Apart from some geometric factors, it describes, through a deficit angle 5^, the parallel 
transport of a vector around an elementary dual lattice loop. It should, therefore, be adequate to 
represent this second term by a sum over contributions over all {d — 3)-dimensional hinges (edges 
in 3+1 dimensions) h attached to the simplex a, giving, therefore, in three dimensions 



Here 5h is the deficit angle at the hinge h, Ih the corresponding edge length, and = \J g{a) the 
volume of the simplex (tetrahedron in three spatial dimensions) labeled by a. Gij (a) is obtained 




(19) 




(20) 




(21) 
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either from Eq. (|19p or from the lattice transcription of Eq. (jlOp 

Gij^kl{(y) = ^g~^^'^{cr)[gik{cr)gji{cr) + gii{a)gjk{a) - gij{a)gki{a)] , (22) 

with the induced metric gij{(j) within a simplex a given in Eq. (jl6p . The combinatorial factor n^^ 
ensures the correct normalization for the curvature term, since the latter has to give the lattice 
version of / ^/g^R = 2 Ylh ^h^h (in three spatial dimensions) when summed over all simplices a. The 
inverse of n^/i counts, therefore, the number of times the same hinge appears in various neighboring 
simplices and consequently depends on the specific choice of underlying lattice structure; for a flat 
lattice of equilateral triangles in two dimensions, n„h = 1/60 The lattice Wheeler-DeWitt equation 
given in Eq. ([2T|) was the main result of a previous paper [1]. 



4 Explicit Setup for the Lattice Wheeler-DeWitt Equation 

In this section, we shall establish our notation and derive the relevant terms in the discrete Wheeler- 
DeWitt equation for a simplex. From now on we shall focus almost exclusively on the case of 2+1 
dimensions. The basic simplex in this case is, of course, a triangle, with vertices and squared edge 
lengths labelled as in Fig. 1. We set t^i = a, — ^02 ~ c. The components of the metric for 
coordinates based at vertex 0, with axes along the 01 and 02 edges, are 

511 = o, 512 = ^(a + c-6), 522 = c. (23) 

The area A of the triangle is given by 

= — [2{ah + hc + ca)-a^ -I? -c^], (24) 
16 

so the supermetric G*-^, according to Eq. is 



AA 

Thus for the triangle we have 



1 -1 -1^ 
-1 1 -1 
-1 -1 1 



G'' = A I -1 1 -1 h (25) 



Si dsj 

and the Wheeler-DeWitt equation is 



^^'^ dsi dsi \da db ~'~ db dc ~^ dc da) ^ ^ 



Figure 1: A triangle with labels. 




Figure 2: Neighbors of a given triangle. The picture illustrates the fact that the Laplacian A(/^) 
appearing in the kinetic term of the lattice Wheeler-DeWitt equation (here in 2+1 dimensions) 

contains edges a, b, c that belong both to the triangle in question, as well as to several neighboring 
triangles (here three of them) with squared edges denoted sequentially hy si = if ... sq = Iq. 
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Figure 3: A small section of a suitable dynamical spatial lattice for quantum gravity in 2+1 
dimensions. 
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where the sum is over the three vertices h of the triangle. 

In the following sections we will be concerned at some point with various discrete, but generally 
regular, triangulations of the two-sphere, such as the tetrahedron, the octahedron and the icosahe- 
dron. These were already studied in some detail in \22\ I23j . A key aspect of the Regge theory is 
the presence of a continuous, local lattice diffeomorphism invariance, whose main aspects in regard 
to their relevance for the 3 + 1 formulation of gravity were already addressed in some detail in 
[1] in the context of the lattice weak field expansion. Here we will add some remarks about how 
this local invariance manifests itself in the 2 + 1 formulation, and, in particular, for the discrete 
triangulations of the sphere studied later on in this paper. Of some relevance is the presence of 
exact zero modes of the gravitational lattice action, refiecting a local lattice diffeomorphism invari- 
ance, present already on a finite lattice. Since the Einstein action is a topological invariant in two 
dimensions, the relevant action in this case has to be a curvature-squared action supplemented by 
a cosmological constant term. Specifically, part of the results in |24^ |22] can be summarized as 
follows. For a given lattice, one finds for the counting of zero modes 

Tetrahedron (Nq = 4) : 2 zero modes 
Octahedron(A''o = 6) : 6 zero modes 
Icosahedron(A'^o = 12) : 18 zero modes . 

(28) 

Thus if the number of zero modes for each regular triangulation of the sphere is denoted by N^.m., 
then the results can be reexpressed as 

N,,m. = 2iVo - 6 , (29) 

which agrees with the expectation that, in the continuum limit, Nq — )• oo, Nz.m./^o should approach 
the constant value d va. d space-time dimensions, the expected number of local parameters for a 
diffeomorphism. Similar estimates were obtained when looking at deformations of a flat lattice 
in various dimensions |22] . The case of near-flat space is obviously the simplest: by moving the 
location of the vertices around in flat space, one can find a different assignment of edge lengths that 
represents the same flat geometry. This then leads to the d- A'^o-parameter family of transformations 
for the edge lengths in flat space. 

In general, lattice diffeomorphisms correspond to local deformations of the edge lengths about 
a vertex, which leave the local geometry physically unchanged, the latter being described by the 

^ Instead of the combinatorial factor rio-h, one could insert a ratio of volumes Vah/Vu (where Vu is the volume per 
hinge 17 and Vch is the amount of that volume in the simplex o), but the above form is simpler. 
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values of local lattice operators corresponding to local volumes, curvatures, etc. The lesson is that 
the correct count of continuum zero modes will, in general, only be recovered asymptotically for 
large triangulations, where Nq is significantly larger than the number of neighbors to a point in 
d dimensions. With these observations in mind, we can now turn to a discussion of the solution 
method for the lattice Wheeler-DeWitt equation in 2 + 1 dimensions. 

One item that needs to be discussed at this point is the proper normalization of various terms 
(kinetic, cosmological and curvature) appearing in the lattice equation of Eq. ()2ip . For the lattice 
gravity action in d dimensions one has generally the following correspondence 

fd'^x^ ^ Y,Va (30) 

where is the volume of a simplex; in two dimensions it is simply the area of a triangle. The 
curvature term involves deficit angles in the discrete case, 

i J d'^x^R ^ Y.Vh6h (31) 

h 

where 6h is the deficit angle at the hinge h, and V/j the associated "volume of the hinge" |12j . 
In four dimensions the latter is the area of a triangle (usually denoted by Ah), whereas in three 
dimensions it is simply given by the length of the edge labeled hy h. In two dimensions V/j = 1. 
In this work we will focus almost exclusively on the case of 2 + 1 dimensions; consequently the 
relevant formulas will be Eqs. ()30p and (13ip for dimension d = 2. 

The continuum Wheeler-DeWitt equation is local, as can be seen from Eq. (|14p . One can 
integrate the Wheeler-DeWitt operator over all space and obtain 



|- (167rG)2 j d^xA{g) + 2X j d^x^- j d^x^R^ ^ = 



(32) 



with the super-Laplacian on metrics defined as 



A(ff) = Gij^kii^) . (33) 

In the discrete case one has one local Wheeler-DeWitt equation for each triangle [see Eqs. (j20p and 
(1211)1. which therefore takes the form 



J -(167rG)2 A(/2) - K ^5i + 2\A^ \ 

y icA ) 



^ = , (34) 



where A(/^) is the lattice version of the super-Laplacian, and we have set for convenience k = 2 ng- h- 
As we shall see below, for a lattice of fixed coordination number, k is a constant and does not depend 
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on the location on the lattice. In the above expression A(/^) is a discretized form of the covariant 
super-Laplacian, acting locally on the space oi s = l'^ variables. From Eqs. (|26p and (j34p one has 
explicitly 

Note that the curvature term involves three deficit angles 6i, associated with the three vertices of a 
triangle. Now, Eq. (|34p applies to a single given triangle, with one equation to be satisfied at each 
triangle on the lattice. One can also construct the total Hamiltonian by simply summing over all 
triangles, which leads to 



I A A A icA J 



(36) 



Summing over all triangles (A) is different from summing over all lattice sites (i), and the above 
equation is equivalent to 



J -(167rG)2 E^(^^) + 2AE^A - Kq ^ = 

I A A i ) 



(37) 



where q is the lattice coordination number and is determined by how the lattice is put together 
(which vertices are neighbors to each other, or, equivalently, by the so-called incidence matrix). 
Here, q is the number of neighboring simplexes that share a given hinge (vertex). For a flat 
triangular lattice q = 6, whereas for a tetrahedron, octahedron, and icosahedron, one has q = 3, 4, 5, 
respectively. For proper normalization in Eq. (I36p one requires 



as well as 



A 

i 

This last correspondence allows one to fix the overall normalization of the curvature term 

K = 2n^h = -, (40) 
Q 

which then determines the relative weight of the local volume and curvature terms. 



5 Choice of Coupling Constants 

As in the Euclidean lattice theory of gravity, we will find it convenient here to factor out an overall 
irrelevant length scale from the problem and set the (unsealed) cosmological constant equal to one 
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as was done in [T7]. Indeed recall that the Euclidean path integral weight always contains a factor 
P{V) oc exp(— AqI^), where V = J ^ is the total volume on the lattice, and Aq is the unsealed 
cosmological constant. The choice Aq = 1 then fixes this overall scale once and for all. Since 
Aq = 2A/167rG, one then has A = SttG in this system of units. In the following we will also find it 
convenient to introduce a scaled coupling A defined as 

A/ 1 ^2 



^"Hi^J '''' 

so that for Aq = 1 (in units of the UV cutoff or, equivalently, in units of the fundamental lattice 
spacing) one has A = l/647rG. One can now rewrite the Wheeler-DeWitt equation so that the 
kinetic term (the term involving the Laplacian) has a unit coefficient and write Eq. (jl4p as 

2A 1 1 

= '''' 

Note that in the extreme strong coupling limit {G — )■ oo) the kinetic term is the dominant one, 
followed by the volume (cosmological constant) term (using the facts about A given above) and, 
finally, by the curvature term. Consequently, at least in a first approximation, the curvature R 
term can be neglected compared to the other two terms in this limit. 

Two further notational simplifications will be done in the following. The first one is introduced 
in order to avoid lots of factors of IGvr in many of the subsequent formulas. Consequently, from 
now on we shall write G 3jS Si shorthand for IGvr G, 

WttG ^G . (43) 

In this notation one then has X = G/2 and A = 1/4G. The above notational choices then lead to a 
much more streamlined representation of the Wheeler-DeWitt equation, 

^ + 7^V9- ^V9r]^ = 0- (44) 



A second notational choice will be dictated later on by the structure of the wave function solutions, 
which will commonly involve factors of For this reason we will now define the new coupling g 

as 

g ^ Vg, (45) 

so that A = 4/(7^ (the latter g should not be confused with the square root of the determinant of 
the metric). 

Later on it will be convenient to define a parameter /3 for the triangulations of the sphere, 
defined as 

P ^ . (46) 

VAG2 
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Factors of 2tt arise here because we are looking at various triangulations of the two-sphere. More 
generally, for a two-dimensional closed manifold with arbitrary topology, one has by the Gauss- 
Bonnet theorem 

I (fx^R = (47) 

with X as the Euler characteristic of the manifold. The latter is related to the genus g (the number 
of handles) via x = 2 — (note that for a discrete manifold in two dimensions one has the equivalent 
form due to Euler x = -^o — -^i + -^21 where Ni denotes the number of simplices of dimension i). 
Thus for a general two-dimensional manifold we will define 

/3 = . (48) 

VAG2 

Equivalently, using 

71 = • = i G^/^ (49) 

and then making use of the coupling g, one has simply 

47r 

H = — (50) 

r 



for the sphere, and in the more general case 



/? = . (51) 



6 Outline of the General Method of Solution 

It should be clear from the previous discussion that in the strong coupling limit (large G) one can, at 
least at first, neglect the curvature term, which can then be included at a later stage. This simplifies 
the problem quite a bit, as it is the curvature term that introduces complicated interactions between 
neighboring simplices (this is evident from the lattice Wheeler-DeWitt equation of Eq. ()2ip . where 
the deficit angles enter the curvature term only). 

The general procedure for finding a solution will be as follows. First a solution will be found for 
equilateral edge lengths s. Later this solution will be extended to determine whether it is consistent 
to higher order in the weak field expansion. Consequently we shall write for the squared edge 
lengths 

E = s{l + ehi^) , (52) 
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with e a small expansion parameter. Therefore, for example, in Eq. ()35p one has a = s(l + eha), 
b = s{l + ehf)) and c = s(l + e/ic). The resulting solution for the wave function will then be 
given by a suitable power series in the h variables. Nevertheless, in some rare cases (such as the 
single-triangle case described below or the single tetrahedron in 3 + 1 dimensions [Ij), one is lucky 
enough to find immediately an exact solution, without having to rely in any way on the weak field 
expansion. 

To lowest order in h, a solution to the Wheeler-DeWitt equation is readily found using the 
standard power series (or Frobenius) method, appropriate for the study of quantum mechanical 
wave equations. In this method one first obtains the correct asymptotic behavior of the solution 
for small and large arguments and later constructs a full solution by writing the remainder as a 
power series or polynomial in the relevant variable. Of some importance in the following is the 
correct determination of the wave functional ^ for small and large areas (small and large s), and 
to what extent the resulting wave function can be expressed in terms of invariants such as areas 
and curvatures, or powers thereof. 

In the following we will see that the natural variable for displaying results is the scaled total 
area x, defined as 

X = iVlAtot = 2\/X Aa . (53) 

A 

We will look at a variety of two-dimensional lattices, including the regular triangulations of the two- 
sphere given by the tetrahedron, octahedron and icosahedron, as well as the case of a triangulated 
torus with coordination number six. In the equilateral case the natural variable for displaying the 
results is then 

X = 2\f\Atot = 2NA\f\AA. (54) 
Later on we will be interested in taking the infinite volume limit, defined in the usual way as 



iVA - 


-> oo , 


Atot - 


oo , 


Atot 




-)■ const 


Na 





(55) 



It follows that this last ratio can be used to define a fundamental lattice spacing Iq, for example 
via Atot/NA = Aa = V3l^/4. 

The full solution of the quantum mechanical problem will, in general, require that the wave 
functions be properly normalized, as in Eq. (jl3p . This will introduce at some stage wave function 
normalization factors M and M, which will be fixed by the standard rules of quantum mechanics. 
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If the wave function depends on the total area only, then the relevant requirement becomes 

/ dAtot I ^(Atot) P = 7^ / dx I ^(x) |2 = 1 . (56) 

-'0 2 V A -'o 

As in nonrelativistic quantum mechanics, two solutions are expected, only one of which will be 
regular as the origin and thus satisfy the wave function normalizability requirement. 

At this point it will be necessary to discuss each lattice separately in some detail. For each 
lattice geometry, we will break down the presentation into four separate items: 

(a) Equilateral Case in the Strong Coupling Limit (e = 0). This subsection will find a 
solution in the extreme strong coupling limit (large G), without curvature term in the Wheeler- 
DeWitt equation. The solution will not rely on the weak field expansion, and the results will be 
exact to zeroth order in the weak field expansion of Eq. (|52|) . In this case the simplices are all taken 
to be equilateral, and the lattice edge lengths fluctuate together. 

(b) Equilateral Case with Curvature Term (e = 0). Next, the curvature term is included. 
The solution again will not rely on the weak field expansion, and all the triangles will be taken to 
be equilateral. The resulting solution will, therefore, be valid again (and exact) to zeroth order in 
the e expansion parameter of Eq. (j52p . 

(c) Large Area in the Strong Coupling Limit (e 7^ 0). In this case we will look at nonzero 
local fluctuations in Eq. (j52p . The method of solution will now rely on the weak field expansion 
for large areas (large s), but nevertheless it will turn out that an exact solution can be found in 
this case. The resulting answer will be found to be correct to arbitrarily large order ©(e"), with n 
a positive integer. 

(d) Small Area in the Strong Coupling Limit (e 7^ 0). Finally we will look at the case of 
nonzero fluctuations [e 7^ in Eq. (|52p ] in the limit of small areas (small s). In this limit we will 
find that, in general, the solution can be written entirely in terms of invariants involving total areas 
and curvatures only up to order 0{e) or ©(e^), depending on whether a further symmetrization of 
the problem is performed or not. 

If the reader is not interested in the details for each lattice, he can skip the next few subsections 
and go directly to the summary presented in Sec. (j6.6p . 

6.1 Single Triangle Case 

From Eq. (j34p the Wheeler-DeWitt equation for a single triangle reads 
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where a, b, c are the three squared edge lengths for the given triangle, and A/^ is the area of the 
same triangle. Note that for a single triangle there can be no curvature term. Equivalently one 
needs to solve 

+ 7^ + 7^ + A^aAc) = 0. (58) 



^ da db db dc dc do 
If one sets 

^[a,&,c] = $[^a] , (59) 
then one finds the following equivalent differential equation 

For a single triangle the total area equals the area of the single triangle, Atot = ^A- Here it will 
be convenient to define 

X = \\f\Atot = 4\/XAa (61) 

so that the solution will be function of this variable only. Note that in this case, and in this case 
only, we will deviate from the general definition of the variable x given in Eq. (j53p . One can then 
write the solution to Eq. ()60p in the form 

^>{x)=N^-^ (62) 



n = - (63) 



Ji/2 ( 4VAA 



Hot 



*(^) = — rirf ■ (64) 



with 



so that 



Hot 

The wave function normalization constant is given here by 

M = 2 A^/^ . (65) 

Note that the above solution is exact, and did not require in any way the weak field expansion. 
Two alternate forms of the wave function are 

sin {^VxAtot 

^{Atot) = M 



2^f2^\J\A, 



Hot 



AAa/- exp(-4i VA^tot) iFi (l,2,8i VAAtot) . (66) 
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Here iFi(a,b, z) is the confluent hyper geometric functions of the first kind. The usefulness of the 
latter representation will become clearer later, when other lattices are considered and the curvature 
term is included. Expanding the solution for small area one obtains 



1 4 
X X 



6 120 ^ ^ 



(67) 



which shows that it is indeed nonsingular and, thus, normalizable. 

In the limit of large areas, a solution to the original differential equation is given either by the 
asymptotic behavior of the above Bessel (here sine) function (J), the same limiting behavior for 
the corresponding Bessel function Y, or by the two corresponding Hankel functions (H). 

^' ~ — exp{±ix) ~ — — exp (±4i\/jAtot) . (68) 

Nevertheless among those four solutions, only one is regular and, therefore, physically acceptable. 

The calculation for a single triangle can be regarded as a useful starting point, and a basic 
stepping stone, for the strong coupling expansion in It shows the physical characteristics 

of the wave function solution deep in the strong coupling regime: for G — )• cxd the coupling term 
between different simplices, which is caused mainly by the curvature term, disappears and one ends 
up with a completely decoupled problem, where the edge lengths in nonadjacent simplices fluctuate 
independently. 



6.2 Tetrahedron 

In the case of the tetrahedron one has 4 triangles, 6 edges, and 4 vertices, and 3 neighboring 
triangles for each vertex. Let us discuss again, here, the various cases individually, 
(a) Equilateral Case in the Strong Coupling Limit (e = 0) 

We first look at the case when e = in Eq. ([52]) , deep in the strong coupling region and without 
the curvature term. 

Following Eq. ([53]) we define the scaled area variable as 

X = 2\f\Atot = 4 X 2\/XAa (69) 

and the solution will be found later to be a function of this variable only. For equilateral triangles 
the wave function ^ needs to satisfy 

2 

+ _ ^' + ^ = . (70) 
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The correct solution can be written in the form 



with 



^>{x)=N^-^ (71) 



n = \ (72) 



Ji/2 {2VlAt. 



so that 

'2\f\Atot 

The wave function normahzation constant is given by 



^(^) = ^ . (73) 



= \/2 A3 . (74) 



Below are two equivalent forms of the same wave function 

sin {2\f\Atot 

*(^tot) = M 



2^^/\Atot 

M \ - e^v{-2iVlAtot) iFi{l,2Ai^\Atot) (75) 



for the equilateral case. In the limit of small area one obtains 

" (76) 



1 4 

1 _ ^ + ^ + o 

6 120 ^ ' 



which again confirms that the wave function is regular at the origin. Since one is solving a second- 
order linear differential equation, one expects two solutions; here, one is singular and the other 
one is not, as is often the case in quantum mechanics. For the geometry of the tetrahedron, one 
solution can be written in terms of Bessel functions of the first kind (J) 

J\l2{x) ^ /2^sinx ^^^^ 

The Bessel function of the second kind (y) also satisfies the same differential equation, but since 

^1/2(3;) _ /2"cosx 



(78) 

'X V TT X 

this second solution is not normalizable, it is therefore discarded on physical grounds. We shall see 
below that the same behavior at small x holds also for the nonzero curvature term. Note that both 
of the above solutions are real. 



^ There are also linear combinations of Bessel functions which give complex Hankel (if) functions. These satisfy 
the Wheeler-DeWitt equation as well; however, they are not physically acceptable since both are singular at the 
origin. 
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(b) Equilateral Case with Curvature Term (e = 0) 

Next we include the effects of the curvature term. To zeroth order in weak field expansions, 
when all edges fluctuate in unison, one now needs to solve the ordinary differential equation 

2 2/3 

*" + - ^ + ^ = , (79) 

X X 

with (3 = In I \f\G'^ as in Eq. (|46p . Since the deficit angle = vr at each vertex, the curvature 
contribution for each triangle is k • vr • 3. In this case one has, therefore, 

_ 1 

and, therefore, the solution is given by 

^ ~ exp(-2iV^Atot) iFi 2, ^i^\At^ 

27r 



2i\r\Ato^ iFi 1-i 2, 4i\/XAtoi (81) 



= exp 

in the equilateral case, up to an overall normalization factor. Note that in this case one had to 
include a factor of Atotj i^'^-^is) (which in the tetrahedron case equals one) in the imaginary part of 
the first argument of \F\. 

(c) Large Area in the Strong Coupling Limit (e / 0) 

Next we look at the case e 7^ in Eq. (j52p . In the limit of large areas one finds that the two 
independent solutions reduce to 



^ ~ exp(±i2;) ~ exp ±2iVA^toi (82) 

to all orders in e. To show this, one sets ^ = e""^*°*, where Atot is a sum of the four triangle areas 
that make up the tetrahedron, and then expands the edge lengths in the usual way according to 
Eq. (|52p . by setting a = s(l + eha), etc. Here we are interested specifically in the limit when s is 
large and e is small. One then finds that the rhs of the lattice Wheeler-DeWitt equation is given 
to 0(e") by 

pQ. "\/3 S "I , / \ 71 

— ^ a" (a2 + 4A) e"s" (V /i) +••• . (83) 

One concludes that in this limit it is sufficient to have 

a2 + 4A = 0, (84) 

ox a = ± li , to obtain an exact solution in the limit n — ^ 00. Note that in the strong coupling 
limit the two independent wave function solutions in Eq. ()82p completely factorize as a product of 
single-triangle contributions. 
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(d) Small Area in the Strong Coupling Limit (e 7^ 0) 

In the limit of small area, we have shown before that the solution reduces to a constant in the 
equilateral case [O (e^)] for small x or small areas. Beyond the equilateral case one can write a 
general ansatz for the wave function in terms of geometric invariants 




and then expand the solution in e for small s. To zeroth order in e we had the solution ^ ~ J„(x)/x"' 
with X = 2 \f\Atot and n = 1/2. This gives in Eq. (j85]) 70 = 0, 72 = — | A and 74 = ^ A^. To 
linear order [O (e)] one finds, though, that terms appear which cannot be expressed in the form 
of Eq. (j85p . But one also finds that, while these terms are nonzero if one uses the Hamiltonian 
density (the Hamiltonian contribution from just a single triangle), if one uses the sum of such 
triangle Hamiltonians, then the resulting solution is symmetrized, and the corrections to Eq. (|85|) 
are found to be of order O (e^) . Then the wave function for small area is of the form 

1 - |AAL + ^A^aL + ... (86) 

up to terms O(e^). 



6.3 Octahedron 

The discussion of the octahedron proceeds in a way that is similar to what was done before for 
the tetrahedron. In the case of the octahedron one has 8 triangles, 12 edges and 6 vertices, with 4 
neighboring triangles per vertex. Again we will now discuss the various cases individually, 
(a) Equilateral Case in the Strong Coupling Limit (e = 0) 

Again we look first at the case e = in Eq. ([52]) , deep in the strong coupling region and without 
the curvature term. Following Eq. (j53|) we define the scaled area variable as 

x = 2\/X^tot = 8x2\/XAa (87) 

and it is found that the solution is a function of this variable only. For equilateral triangles the 
wave function ^ needs to satisfy 

4 

+ + $ = . (88) 

X 

The correct solution can be written in the form 

^(-) = ^ (89) 
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with 



so that 



n 



J3/2 {2 VxAtot 



2 VlAtot) 

The wave function normahzation factor is given by 



3/2 



(90) 



(91) 



Equivalent forms of the above wave function are 



T5AV4. 



(92) 



^!{A 



tot 



exp (^-2i\/lAtot) iFi (^2,4,4 i\/A^tot) 
sin (2 71 A 



23/2 p (5) 

COS (2 \f\Atot 



2J2^\Al 



+ 



tot 



4V2^A3/2 Al 



tot 



(93) 



These can be expanded for small At^t or smah x to give 



^ = M 



V2 
2,^ 



2 4 
10 280 ^ ' 



(94) 



We note here again that both Bessel functions of the first (J) and second (Y) kind, in principle, 
give solutions for this well as the two corresponding Hankel (H) functions. Nevertheless, 

only the solution associated with the Bessel J function is regular near the origin, 
(b) Equilateral Case with Curvature Term (e = 0) 

Next, we include the effects of the curvature term. Since here the deficit angle 6 = 2tt/3 at 
each vertex, the curvature contribution for each equilateral triangle is k • ^ • 3 = 27r k. For the 



octahedron one has in Eq. (j40 



l^octa — 2 • 



(95) 



With the curvature term one finds 

*(Atot) ^ exp 2 i VI Atot J i-ri 

(^-2i VI Atot) iFih-i 



F, 4, 4iVlA, 



exp 



VI 
2tt 



tot 



4, AiV^A 



Hot 



(96) 



Note that in this case one had to include a factor Atot/i^A/^), which in the octahedron case equals 



two. 



(c) Large Area in the Strong Coupling Limit (e / 0) 
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In the limit of large areas the two independent solutions reduce to 

~ exp{±ix) ~ exp (±2iVJAtot) (97) 

to all orders in e. In other words, to O (e") with n — )• oo, as for the tetrahedron case. Note also 
that in the strong coupling limit the two independent wave function solutions again completely 
factorize as a product of single-triangle contributions, 
(d) Small Area in the Strong Coupling Limit (e 7^ 0) 

In the limit of small area, the solution approaches a constant in the equilateral case. To go 
beyond the equilateral case, one can write again a general ansatz for the wave function, written in 
terms of geometric invariants as in Eq. (|85|) . Then the solution can be expanded in e for small s. 
To zeroth order in e, the solution is ~ Jn{x)/x^ with n = 3/2. This gives in Eq. (j85p 70 = 0, 
72 = — I A and 74 = ^ A^. However, to linear order [O (e)] one finds again that linear terms in 
h appear which cannot be expressed in the form of Eq. (j85p . But one also finds that while these 
terms are nonzero if one uses the Hamiltonian density (the Hamiltonian contribution from just 
a single triangle), if one uses the sum of such triangle Hamiltonians then the resulting solution 
is symmetrized, and the corrections to Eq. (|85p are found to be of order O (e^) . Then the wave 
function for small area is of the form 

^c.l-fXAl, + l^~X'Ai, + ... (98) 

up to terms of 0{e). 
6.4 Icosahedron 

The discussion of the icosahedron proceeds in a way that is similar to what was done before for the 
other regular triangulations. Here one has 20 triangles, 30 edges and 12 vertices, with 5 neighboring 
triangles per vertex. Let us again discuss the various cases individually, 
(a) Equilateral Case in the Strong Coupling Limit (e = 0) 

Again we look first at the case e = in Eq. (|52p , deep in the strong coupling region and without 
curvature term. Following Eq. (|53p we define the scaled area variable as 

X = 2 \/X Atot = 20x2 (99) 

and a solution is found which is a function of this variable only. For equilateral triangles the wave 
function ^ needs to satisfy 

5'" + — 1-' + I- = . (100) 

X 
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A solution can then be found of the form 



with 



so that 



=M'^ (101) 



n = - (102) 



J9/2 I 2 y \Atot 

^(^) = . (103) 



(2VlAtot) 

The wave function normahzation factor is given by 

M = 9 V12T55 A^/^ . (104) 
Below is an equivalent form of the same solution 

^{Atot) = M ^3^^ ^ ^11^ exp [-2iV\Atot) iFi (5, 10, iiVlAtot) . (105) 



For small area Atot or small x, one obtains 

1 



^ = M 



2 4 

^-i2+i^+^(^^) 



(106) 



29/2 r (U) 

which shows that the above solution is regular at the origin and normalizable. 

(b) Equilateral Case with Curvature Term (e = 0) 

Next we include again the effects of the curvature term. Since now the deficit angle J = 7r/3 at 
each vertex, the curvature contribution for each triangle is k • ^ • 3 = vr k. For the icosahedron one 
has in Eq. (l40l) 

l^icosa — 2 • — . (10''') 
5 

Then with the curvature term included for equilateral triangles one obtains for equilateral triangles 

Pie')] 

^{Atot) ^ e^p(^-2iVJAtot) iFi(^5-i^^^^^, W, A^^/lAto^ 

= e^p 2 i VI Atot) iFi (^5- i-^^, 10, Ai VI Atot^ ^ (108) 

up to an overall wave function normalization constant. Note that in this case one had to include a 
factor Atot/^A^, which in the dodecahedron case equals five. 

(c) Large Area in the Strong Coupling Limit (e 7^ 0) 

26 



In the limit of large areas the two independent solutions reduce to 

Atot (109) 

to all orders in the weak field expansion parameter e, as for the tetrahedron and octahedron case. 
Note also that in the strong coupling limit the two independent wave function solutions again 
completely factorize as a product of single-triangle contributions, 
(d) Small Area in the Strong Coupling Limit (e 7^ 0) 

In the limit of small area, the solution approaches a constant in the equilateral case. To go 
beyond the equilateral case, one can write again a general ansatz for the wave function, written in 
terms of geometric invariants as in Eq. (|85p . Then the solution in e for small s. To zeroth order 
in e the solution is ^' ~ Jn{x)/x^ with n = 9/2. This gives in Eq. ([85]) 70 = 0, 72 = — ^ A and 
74 = A^. But to linear order [O (e)] one finds again that linear terms in h appear which cannot 
be expressed in the form of Eq. (jSSp . But one also finds that, while these terms are nonzero if 
one uses the Hamiltonian density (the Hamiltonian contribution from just a single triangle), if one 
uses the sum of such triangle Hamiltonians then the resulting solution is symmetrized, and the 
corrections to Eq. (|85p are found to be of order O (e^) . Then the wave function for small area is of 
the form 

^ ^ 1 - ^aaL + m>^^<t + ■■■ , (110) 

up to terms of 0{e). 



6.5 Torus 

Finally we will consider a regularly triangulated torus, which will consist here of an infinite lattice 
built out of triangles, with each triangle having 12 neighboring triangles. The torus topology is 
equivalent to requiring periodic boundary conditions in the two spatial directions. Of course, one 
could consider the same type of lattice but with some other sort of boundary condition, but we 
shall not pursue that aspect here. 

Due to the local structure of the lattice Wheeler-DeWitt equation in Eq. ()34p . it will not be 
necessary to include in the wave function triangles that are arbitrarily far apart. Instead it will be 
sufficient, in order to determine the overall structure of the solution, to include only those triangles 
that are affected in a nontrivial way by the interaction terms in the Wheeler-DeWitt equation. In 
the present case, this requires the consideration of one given triangle plus its 12 neighbors, giving 
a total of 13 triangles. Here, we will also set as before x = 2 y^Atot- 
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(a) Equilateral Case in the Strong Coupling Limit (e = 0) 

For this case the relevant equation and its solution are largely in line with what was obtained 
for the previous cases. For equilateral triangles the wave function ^ has to satisfy 

13 

V + — + * = . (Ill) 
The wave function can now be written as 

^{x)=M^ (112) 



X" 



with, here, (due to our specific choice of sublattice) 



4 

so that 

•^11/4 



n = H (113) 



(2VxAtot) 

The wave function normalization constant is given in this case by 



= ^f — ^ TTTjf ■ (114) 



For the above wave function an equivalent form is 



^(Atot) = -^^[TJIY^ c^p[- 2 i VI Atot) iFi (^^,^,4iVlAto?j . (116) 



15 570 ^ ^ 



(117) 



Expanding the above solution for small area one obtains 

^ = ■'^ 211/4 r (^) 

which shows the above solution is indeed regular at the origin. 

(b) Equilateral Case with Curvature Term (e = 0) 

In the case of the torus, the curvature term is zero (x = 0), so there are no changes to the 
preceding discussion. 

(c) Large Area in the Strong Coupling Limit (e ^ 0) 

In the limit of large areas the two independent solutions reduce to 



a; — > oo 



exp(±za;) ~ exp (^±i2\/jAtot^ (118) 



to all orders in e. This is similar to what was found earlier for the other lattices. In particular, the 
two independent solutions again completely factorize as a product of single-triangle contributions. 
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(d) Small Area in the Strong Coupling Limit (e 7^ 0) 

In the limit of small area, the regular solution approaches a constant and the discussion, and 
solution, is rather similar to the previous cases. Here one finds 

* :^ 1 - ^A^L + Js^'^t + ••• ' (119) 

up to terms of O(e^). 

6.6 Summary of Results 

In this section we will summarize the results obtained so far for the various finite lattices considered 
(tetrahedron, octahedron, icosahedron, and regularly triangulated torus), 
(a) Equilateral Case in the Strong Coupling Limit (e = 0) 

It is rather remarkable that all of the previous cases (except the trivial case of a single triangle, 
which has no curvature) can be described by one single set of interpolating wave functions, where 
the interpolating variable is simply related to the overall lattice size (specifically, the number of 
triangles) . 

Indeed for equilateral triangles and in the absence of curvature, the wave function "^(x) for all 
previous cases is a solution to the following equation 

2r? -I- 1 

^// ^ ^r' + * = , (120) 

X 

with parameter n given by 

n = 1{Na-2) (121) 
where N\ = N2 is the total number of triangles on the lattice. Thus 

and, consequently, 

n-tetrahedron 
'^octahedron 
f^icosahedron 
f^torus 

Note that for a single triangle one has n = ^ as well, but the definition of the scaled area is different 
in that case. 



= 1(4-2) = 5, 

3 

2 ' 



2) 



= i(20-2) = - , 
= i(13-2) = - 



(122) 



(123) 
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Furthermore the differential equation in Eq. ()120p describes, in spherical coordinates and with 
suitable choice of constants, the radial wave function for a free quantum particle in D = 2n + 2 
dimensions. Indeed recall that in D dimensions the Laplace operator in spherical coordinates has 
the form 

d"^^ D-l I ^ ^ , 

Al- = — ^ + -jr + - ^s^-' * (124) 

ar^ r or 

where Ago-i is the Laplace-Beltrami operator on the (D — l)-sphere. In our case, the wave function 
does not, to this order, depend on angles and therefore the last (angular variable) term does not 
contribute. The role of the angles is played in our case by the h variables, which to this order do 
not fluctuate. 

A nonsingular, normalizable solution to Eq. (jl20p is then given by 

^{x) = M^^^ = AAe"*^ iFi (n + i,2n + l,2zx) (125) 
where M is the wave function normalization constant 

1/2 



and 



A^/^ , (126) 



^ = IT-^ ttAA. (127) 

2"T(n + l) ^ ^ 

Here and in Eq. (jl25p i-Pi(a, b; z) denotes the confluent hypergeometric function of the first kind, 
sometimes denoted also by M(a, b,; z). In either form, the above wave function is real, in spite of 
appearances. The general asymptotic behavior of the solution ^(x) is found from Eq. (jl20p . For 
small X one has 

^(x) ~ x" (128) 

with index a = 0, —2n. The latter solution is singular and will be discarded. For large x one finds 
immediately 

^(x) exp(±ix), (129) 

X 2 

which is of course consistent with all the previous results. Indeed the other possible independent 
solution of Eq. (jl20p would be 

^(x) ^ , (130) 

where Yn{x) is a Bessel function of the second kind (or Neumann function). However, the latter 
leads to a wave function ^ which is singular as x — )• 0, 

^(x) -r(n)2"x-2" (131) 

vr 

30 



and gives, therefore, a solution which is not normahzable. For completeness we record here the 
small X (small area) behavior of the normalized wave function in Eq. (jl25p 

^(x) ~ M ——^ -, (132) 

^ ^ 2'^r(n + l)' ^ ^ 

and the corresponding large x (large area) behavior 

^(x) ~ AAy|-l^cos(x-f -9 , (133) 

X 2 

both of which reflect well-known properties of the Bessel functions Jn(x). 
(b) Equilateral Case with Curvature Term (e = 0) 

When the curvature term is included in the Wheeler-DeWitt equation, and still in the limit of 
equilateral triangles, one obtains the following interpolating differential equation 

„ 2n + 1 , 2/3 , , 

^" + ¥ ^^r + ^r = 0, 134 

X X 

which now describes the radial wave function for a quantum particle in D = 2 n + 2 dimensions, 
with a repulsive Coulomb potential proportional to 2/3. The nonsingular, normalizable solution is 
now given by 

^{x) ~ e""' 1^1 (n + I -i/3, 2n + 1, 2i3;) , (135) 

up to an overall wave function normalization constant M{n,f3). The normalization constant can 

be evaluated analytically but has a rather unwieldy form and will not be recorded here. Note 

that the imaginary part (/3) of the first argument in the confluent hyper geometric function of 

Eq. ()135p depends on the topology but does not depend on the number of triangles. In view of 

the previous discussion the parameter n increases gradually as more triangles are included in the 

simplicial geometry. For the regular triangulations of the sphere, the total deficit angle (the sum 

of the deficit angles in a given simplicial geometry) is always in, so even if one writes for the wave 

functional "^lAtot, Stot], the curvature contribution J2h^h ^ constant and does not contribute in 

any significant way. Note also that, in spite of appearances, the above wave function is still real for 

nonzero f3. That ^(x) in Eq. (jl35p is a real function can be seen, for example, from its definition 

via the power series expansion 

2/3 l + 2n-4^2 „ /3 (5 + 6n - 4/32) „ . ^ ^ 

^(x) ~ 1 + —X ^x^ -— ^ ^ ' x^ + O(x^) , (136) 

^' 2n+l 4+12n + 8n2 6 (3 + lln + 12n2 + 4^3) ^ ^' ^ ^ 

and again up to an overall normalization factor J\f{n,f3). 

The general asymptotic behavior of the solution ^{x) is again easily determined from Eq. (jl34p . 

For small x one has 

^(x) ~ x" (137) 
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with again a = 0, —2n, and, therefore, independent of the curvature contribution involving /3. The 
second solution is singular and will be discarded as before. For large x one finds immediately 

^{x) ^- exp{±i(x-/31nx)}, (138) 

which is of course consistent with all previous results. It also shows that the convergence properties 
of the wave function at large x are not affected by the /3 term. A second independent solution to 
Eq. ()134p is given by 

^(x) ~ e''"" U (n + i-i/3, 2n + l, 2ix) , (139) 

where U (a, b, ; z) is the confluent hypergeometric function of the second kind (sometimes referred 
to as Tricomi's function). This second solution is singular at the origin, leading to a wave function 
that is not normalizable and will not be considered further here. 

The asymptotic behavior of the regular solution for large argument z (discussed in standard 
quantum mechanics textbooks such as \25\ [26] and whose notation we will follow here) can be 
obtained from the asymptotic form of the confluent hypergeometric function iFi, defined originally, 
for small z, by the series 

,-, / , N -, a(a + l)z^ ^. 

It is common procedure to then write i-Fi(a, b, z) = Wi{a, b, z) + W2{a, 6, z), where Wi and W2 are 
separately solutions of the confluent hypergeometric equation 

d'^F dF 

Then an asymptotic expansion for iFi (or M) is obtained from the following relations: 

Wi{a,b,z) = -S^{-zr-w{a,a-b+l,-z) (142) 
L [b — a) 



where 



W2(a,b,z) = — f^e^z" w(l - a,b - a,z) 
r a 



w{a,(3,z) ~ 1 + ^ + ^ ^ '- + ■■■, 143 

z -> 00 2; 1! 2! 



with the irregular (at the origin) solution given instead by the combination G{a, b, z) = i Wi (a, b, z) — 
iW2{a,b, z). One immediate and useful consequence of the above result is that, as anticipated be- 
fore, the behavior of the regular solution close to the origin is not affected by the presence of the 
j3 (curvature) term. In other words, the wave function solution ^{x) in Eq. ()135p is always well 
behaved for small areas and, therefore, leads to a perfectly acceptable, normalizable solution. 
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Furthermore, the combination and properties of arguments in the confluent hypergeometric 
function in Eq. (|135p allow one to write it equivalently as a Coulomb wave function with (Sommer- 
feld) parameter rj 



Ci{r,) p'+^ -e-'PiFiil + l- ir],2l + 2,2ip) = Fiir],p), 



(144) 



where Fi{r], p) denotes the regular Coulomb wave function that arises in the solution of the quantum 
mechanical three-dimensional Coulomb problem in spherical coordinates |25^ [26] . The latter is a 
solution of the radial differential equation 



dp^ 



+ 



\_2ri _ + 



Fi = 



(145) 



p p- 

with the actual radial wave function then given by Ri{r) = Fi{kr)/r. After comparing the above 
equation with Eq. (jlSSp one then identifies p = x, / = n — ^ and rj = j3. Thus I = — 1, where 

A^A is the number of triangles on the lattice. The proportionality constant Ci in Eq. ()144p is given 
by the (Gamow) parameter 

2^6"^ |r(/ + 1 + ir])\ 



Ciiv) 



(146) 



r(2/ + 2) 

One then has immediately, from Eq. ()135p . an equivalent representation for the regular wave func- 
tion as 



1 



(147) 



again up to an overall wave function normalization constant J\f{n,f3). Again we note here that, 
on the other hand, the irregular Coulomb wave function [usually denoted by Gi{rj,p)] is singular 
for small r and will, therefore, not be considered here. Further relevant properties of the Coulomb 
wave function can be found in [251 12511^ . 

The known asymptotics of Coulomb wave function |27i I29j allow one to derive the following 
result for the wave function ^ for large x 

1 



with (Coulomb) phase shift 
Also, from Eq. (fTi6D . 



sm 



C„ i(/3)-x' 

n 2 



n+ 



(2n - l)7r 
a; — p In 2x h an 



(148) 



UniP) = argr(n + i+i/3) . 



2"-te-'^|r(n + l + 
T(2n + 1) 



C 



2 



(149) 



(150) 
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It is easy to check that the above result correctly reduces to the asymptotic expression given earlier 
for \I' in Eq. (jl33p in the limit /? = 0. The structure of the wave function in Eq. (jl48p implies 
that the norm is still finite for /3 7^ 0, since the convergence properties of the wave function are not 
affected by the curvature term. 

(c) Large Area in the Strong Coupling Limit (e 7^ 0) 

In the limit of large areas the two independent solutions reduce to 

^ ~ exp{±ix) (151) 
x — > 00 

where x oc Atot- This is true without assuming the weak field expansion, as was already the case 
before (see, in particular, the section discussing the tetrahedron case). 

Consequently in the strong coupling limit the two wave function solutions in Eq. p5ip com- 
pletely factorize as a product of single-triangle contributions, 

^' ~ JJ exp (^±2i\/XAA) , (152) 
A 

again up to an overall normalization constant. The above result, anticipated in [T], was the basis for 
the variational treatment using correlated product wave functions given in our previous work. Note 
also, in view of the result of Eq. (jl33p . that the correct solution, satisfying the required regularity 
condition for small areas, is actually a linear combination of the above factorized solutions. 

(d) Small Area in the Strong Coupling Limit (e ^ 0) 

In the limit of small area, we have shown before in all cases that the solution reduces to a 
constant in the equilateral case [O (e'^)] for small x or small areas. To linear order [O (e)] the 
general result is still that linear terms in h appear which cannot be expressed in the form of 
Eq. (|85p . But one also finds that, while these terms are nonzero if one uses the Hamiltonian 
density (the Hamiltonian contribution from just a single triangle), if one uses the sum of such 
triangle Hamiltonians then the resulting solution is symmetrized, and the corrections to Eq. (j85p 
are found to be of order O (e^) . In other words, it seems that some residual lattice artifacts that 
survive at very short distances can be partially removed by a suitable coarse-graining procedure on 
the Hamiltonian density. 

One might wonder what lattices correspond to values of n greater that 9/2, which is the highest 
value attained for a regular triangulation of the sphere, corresponding to the icosahedron. For 
each of the three regular triangulations with Nq sites one has for the number of edges A'^i = ^Nq 
and for the number of triangles A'2 = (f ~ l)-^0 + 2, where q is the number of edges meeting 
at a vertex (the local coordination number). In the three cases examined before q was between 
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three and five, with six corresponding to the regularly triangulated torus. Note that for a sphere 
A'^o — A'^i + = 2 always. The interpretation of other, even noninteger, values of q is then clear. 
Additional triangulations of the sphere can be constructed by considering irregular triangulations, 
where now the parameter q is interpreted as an average coordination number. Of course the 
simplest example is a semiregular lattice with Na vertices with coordination number qa and Nh 
vertices with coordination number qi,, such that Na + Nb = Nq. Various irregular and random 
lattices were considered in detail some time ago in [16], and we refer the reader to this work for a 
clear exposition of the properties of these lattices. 

We conclude this section by briefly summarizing the key properties of the gravitational wave 
function given in Eqs. ()135p and (]147p . which from now on will be used as the basis for additional 
calculations. First we note that the above wave function is a function of the total area and total 
curvature only and, as such, is manifestly diffeomorphism-invariant and in accord with the spatial 
difFeomorphism constraint. While it was derived by looking at the discrete triangulations of the 
sphere, it contains a parameter n, related to the total number of triangles on the lattice by Eq. (jl2ip . 
that will allow us to go beyond the case of a finite lattice and investigate the physically meaningful, 
and presumably universal, infinite volume limit n — t- oo [see Eq. (j55|) ]. We have also shown that the 
above wave function is, in all cases, an exact solution of the full lattice Wheeler-DeWitt equation of 
Eq. ()2ip in the limit of large areas, and to all orders in the weak field expansion. Again, this last case 
is most relevant for taking the infinite volume limit, defined previously in Eq. (|55p. Furthermore, 
the small area behavior of the wave function plays a crucial role in uniquely constraining, through 
the regularity condition, the correct choice of solution. In this last limit one also finds that the 
various individual lattice solutions agree with the universal form of Eqs. (|135p and (|147p only 
to a low order in the weak field expansion, which is expected given the different short distance 
lattice artifacts of the regular triangulation solutions. Nevertheless, knowledge of their behavior 
is completely adequate for extracting the most important physically relevant piece of information, 
namely the constraint on the wave function based on the stated regularity condition at small areas, 
which comes down to a simple integrability or power counting argument. 

7 Average Area 

In this section we will look at a natural quantum mechanical expectation value, the average total 
physical area of the lattice simplicial geometry. It is one of many quantities that can be calcu- 
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Figure 4: Wave Function versus total area for the octahedron lattice, with and without curvature 
contribution. The wave function is shown here for g = \fG = 1, a value chosen here for illustration 
purposes. The relevant expression for the wave function is given in Eq. ()153p . We refer to the text 
for further details on how the wave function was obtained, and what its domain of validity is. The 
wave functions shown here have been properly normalized. Note that with a nonzero curvature 
term the peak in the wave function moves away from the origin. 

lated within the lattice quantum gravity formalism, and is clearly both manifestly geometric and 
diffeomorphism-invariant. Here we will use the wave functions given in Eqs. ()135p and (|147p . origi- 
nally obtained for the tetrahedron, octahedron and icosahedron, and later extended to any number 
of triangles A'^a 

^{Atot) ^ e-*^iFi (^n + i-i/3,2n + l,2i^^ , (153) 

with n = J {Na — 2), /? = ATr/g^ and g = and again valid up to an overall wave function 
normalization constant. Due to the structure of the wave function the resulting probability dis- 
tribution for the area is rather nontrivial, having many peaks associated with the infinitely many 
minima and maxima of the hyper geometric function. Clearly the most interesting limit is one 
where one considers an infinite number of triangles, — )• oo, which corresponds to n — )• oo in 
Eq. ()153p . In Figs. [l]and[5l we display the behavior of the wave function in Eq. (jl53p . both with 
and without the curvature contribution in the Wheeler-DeWitt equation. One notices that when 
the curvature term is included (/3 7^ 0), the peak in the wave function shifts away from the origin. 
This is largely expected, based on the contribution from the repulsive Coulomb term in the wave 
equation of Eq. (jl34p . 
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Figure 5: Same wave function ^ as in Fig. HI but now for the icosahedron lattice. 

The average total area can then be computed from the above wave function, as the ground state 
expectation value 

(^Ul^) _ I df,[g]A{g)\^{g)\^ 

where g, here, is the three-metric, and dfj,[g] denotes a functional integration over all three- metrics. 
In our case we use the measure 

//•oo 
dti[g] dAtot, (155) 

which then gives, in terms of the scaled area variable x, 

dxx-\^{x)\^ , , 

Jo dx\^{x)Y 

In the absence of a curvature term in the Wheeler-DeWitt equation (/3 = 0), the average area can 
easily be computed analytically in terms of Bessel function integrals, and the result is 

^ 7r(4n-l)r(4n-2) 

Note that the average area diverges as n — i, which corresponds to the tetrahedron; this entirely 
spurious divergence prevents us from using the tetrahedron lattice in plotting and numerically 
extrapolating the remaining two lattices (octahedron and icosahedron) to the infinite lattice limit. 
For the octahedron one finds < Atot >= 15 ^/tt, for the icosahedron < Atot >= 21879 5/39207r, and 



in the large n limit < Atot >= \/2n/7r g + 0{l/^/n). 
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One finds that in tlie presence of a curvature term (/3 7^ 0) tlie resulting integrals are significantly 
more complicated. We have, therefore, resorted to a number of tools, which include an analytic 
expansion in /3, the use of known asymptotic expansions for the wave function at large arguments, 
and an exact numerical integration of the resulting integrals. Let us first discuss here the expansion 
in /3. It is known that the Coulomb wave functions can be expanded in terms of spherical Bessel 
functions (Neumann expansion) [271 [2H1 ISH] , so that one has 

2p 



IT 



r(/ + |)Q(r/)p 



00 

E 

k=l 



bk{v)J^^ lip) 



(158) 



with coefficients bkir]) given by a simple recursion relation. When written out explicitly, the ex- 
pression in curly brackets involves 



J. 



2/_+3 

T+T 



v-J, 



1+ f 



3 (x) + ^-r V ■ J. 



l + l 



1+ 



5(2:) + 



(159) 



with additional terms linear in r] reappearing at higher orders. That the above expansion is a bit 
problematic is not entirely surprising, given the modified asymptotic behavior of the Coulomb wave 
functions for 7] ^ 0. In the following, in order to provide initially some insight into the effects of 
the T] (or (3) term on the wave function ^, we will include the first correction as a perturbation, 
and drop the rest. Later on, higher order corrections can be included as additional contributions. 
With this truncation, the Coulomb wave function in Eq. (|144p becomes 



T r ^ 2/ + 3 ^ , , 



(160) 



with the last term treated as a perturbation, giving for the wave function itself [see Eq. ()135p ] 



^{x) ~ e-*^ iFi (n + i - i/3, 2n + 1, 2ix) 

2n + 2 

J„(x) + /3 ^ J„+i(x) + •• 



(161) 



again up to an overall wave function normalization constant M. Note that if m Bessel function 
terms are kept in Eq. (jl6ip . beyond the zeroth order, strong coupling, term involving Jn(x), then 
the resulting expansion in /3 contains terms up to /S™. One finds to lowest order (m = 1) 
1 _ r(n) 4i-'^(n + l)/3 



2r (n + i)r(2n + i) 



+ 



(162) 



(2n+ l)r(n+ 1)2 

From the above expressions, the average area can then be computed as some still rather complicated 
function. 



^(4n-l)r(4n-2) , 4(n + l)/3 

< ^tot >— 9 < „Q„ K^/ ^. H 



28n-5r(n)^ 



2n + 1 



1 



4i-2nr (n_i)r(n + i)r(2n+l) 



6 



+ 



(163) 
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Figure 6: Average area of a single triangle vs. g = \fG for the octahedron and the icosahedron 
configurations. The average area was calculated using the expression in Eq. (1156(1 . Note the 
qualitative change when one includes the curvature term, with a minimum appearing at 5 ~ 

Additional terms can later be included in the Bessel function expansion of Eq. ()158p , so as to obtain 
more accurate values for the averages; this will be done later. 

Figure [6] shows the exact value of the average area for a single triangle < A a >=< Atot > /^^A 
as a function of the coupling g, obtained by doing the integral in Eq. (jl56p numerically, with the 
wave function given in Eq. ()153p . One noteworthy aspect is that a qualitative change seems to 
occur when one includes the curvature term: a well defined minimum occurs at g ^ 1, which would 
suggest the appearance of some sort of phase transition. Doing the integrals numerically one finds 
a minimum in the average area of a triangle at gc ^ 3.1 for the octahedron, and at gc ~ 2.6 for the 
icosahedron. On the other hand, using the lowest order Bessel function expansion of Eq. (I16ip for 
the octahedron (n = 3/2) one finds a minimum at gc = 2.683, and for the icosahedron (n = 9/2) at 
gc = 2.271. Adding one more Bessel function correction term then gives gc = 3.135 and gc = 2.637 
for the two cases, respectively, which suggests that the expansion is converging. 

The limit of a large number of triangles A'^a — ^ 00 corresponds to taking the parameter n in 
Eq. (|153p to infinity, since n = j {Nj\ — 2). From the lowest order Bessel function expansion one 
obtains the following analytic expression for the average total area 

<.4„>=,.71^[l + -f + 0(4)1 +?<^,5 + ..., (164, 

V'7r[ Ion \™/J ^ 

with /3 = 4:7: /g^ [see Eq. (|50p ]. In this limit the resulting function of g has, again, a well defined 
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minimum at 

or ~ for large n with one Bessel function correction term. With two Bessel function 

correction terms in Eq. ()16ip one finds gc — which again suggests that the expansion 

is slowly converging. Using the exact wave function to do the integrals numerically one finds 
for the minimum — which is close to the above answer. Interestingly enough, the 

above result would suggest that in the limit of infinitely many lattice points the critical point gc 
actually moves to the origin, indicating a phase transition located at exactly g = (G = 0) in the 
infinite volume (n — )• oo) limit (see further discussion later). We note here that the average area 
for a single triangle is obtained by simply dividing the average total area by the total number of 
triangles A'^a = 4n + 2, which then gives in the same limit of large n and strong coupling 

<Aa>= £ + 0{-) . (166) 

Quite generally, the average of the area per site in the lattice theory (the spatial volume per site) 
appears to be well defined mainly due to our wave function normalization choices and, consequently, 
can be explicitly calculated without any leftover ambiguity. 

As will be discussed further below in more detail, the estimate for the critical point given in 
Eq. (jl65p is also in good agreement with a previous variational estimate. In [1] the quantum- 
mechanical variational (Rayleigh-Ritz) method was used to find an approximation for the ground 
state wave function, using as variational wave function a correlated (Jastrow-Slater) product of 
single-triangle wave functions. There it was found, from the roots of the equation < > = 0, 

that the variational parameters are almost purely imaginary for strong coupling (large G > Gc), 
whereas for weak enough coupling (small G < Gc) they become real. This abrupt change in 
behavior of the wave function at Gc then suggested the presence of a phase transition. With the 
notation used in this paper, the result of [1] reads gc ~ in qualitative agreement with the 

result of Eq. ()165p . in the sense that both calculations point to a critical point Gc = in the infinite 
volume limit. 

Let us now make some additional comments which should help clarify the interpretation of the 
previous results. It is well known that if there is some sort of continuous phase transition in the 
lattice theory, the latter is generally associated with a divergent correlation length in the vicinity 
of the critical point. In our case it is clear that at strong coupling (large g) the correlation length is 
small (of order one) in units of lattice spacing. This can be seen from the fact that (a) the coupling 
term in the Wheeler-DeWitt equation is due mainly to the curvature term, which is small for large 
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g, and (b) that the ground state wave function is of the form of a correlated product in the same 
Umit [see Eq. (|152p ]. Then as the effects of the curvature term are included, the correlation length 
starts to grow due to the additional coupling between edge variables. The previous calculation 
would then suggest that the point of divergence is located at g = 0. It is, of course, essential that 
one looks at the limit of infinitely many triangles, A^a oo, since no continuous phase transition 
can occur in a system with a finite number of degrees of freedom. 

It is also of interest here to discuss how the above (Lorentzian) results relate to what is known 
about the corresponding Euclidean lattice theory in three dimensions, which was studied in some 
detail in p4]. There a phase transition was found between two phases, with the weak coupling phase 
G < Gc exhibiting a sort of pathological behavior, whereby the lattice collapses into what geomet- 
rically could be described as a branched polymer. This is clearly a nonperturbative phenomenon 
that cannot be seen from perturbation theory in G. In the Euclidean formulation, average volumes 
are obtained as suitable derivatives of log Ziatt with respect to the bare cosmological constant Aq, 
where Ziatt is the lattice path integral 

Ziau = y"[d/V'""^''^ (167) 
with, in four dimensions, the action given by 

Iiatt = ^oYl ^hih - ^ E ^h{lMh{l^) (168) 

h h 

and h denoting a hinge [more details can be found in |24j]. Similarly, the average curvature can 
also be obtained as a derivative of logZ with respect to A; = l/(87rG). More importantly, a 
nonanalyticity in Z, as induced by a phase transition, is expected to show up in local averages as 
well. From the above expression for Ziatt, exact sum rules can be derived relating various averages 
[30j . In the case of the three-dimensional Euclidean theory the sum rule reads 

2\o<Y,yT> -k<Y,^hlh> -Co = Q (169) 

T h 

where the first term contains a sum over all lattice tetrahedra, and the second term involves a sum 
over all lattice hinges (just edges in this case). The quantity Co, here, is a constant that solely 
depends on how the lattice is put together (i.e. on the local coordination number, or incidence 
matrix) . 

In [23] it was found that the average curvature goes to zero at some gc with a characteristic 
universal exponent 5, 

<Y,Shlh > = -Ro \g-gcf (170) 

h 
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and that the curvature fluctuation diverges in the same limit. From the sum rule in Eq. ()169p one 
then deduces that the average volume in the Euclidean theory has a singularity of the type 

<Y^Vt>= Vo-Vi \g-gcf (171) 

T 

with the same exponent 6 ~ 0.77. The latter is related by standard universality and scaling 
arguments [311 [32l [33j (see [19j for details specific to the gravity case) to the correlation length 
exponent ly hy ly = {1 + 6)/d in d dimensions. To compare to the Lorentzian theory discussed 
in this paper, one notes that the three-dimensional Euclidean theory corresponds to the (2 + 1)- 
dimensional Wheeler-DeWitt theory, so that the average volume in the above discussion should be 
taken to correspond to an average area in our case. Ifl To conclude, the results for the average area 
suggest the existence of a phase transition in the Lorentzian theory located at g = 0. In the next 
sections we will present a further test of this hypothesis, based on physical observables that can 
establish directly and unambiguously the location of the phase transition point. 



8 Area Fluctuation, Fixed Point and Critical Exponent 

Another quantity that can be obtained readily from the wave function ^ is the fluctuation in the 
total area 

XA = ^{< (Atotf > - < Atot >'} . (172) 
The latter is related to the fluctuations in the individual triangles by 

XA = Na {<Al> - <Aa>^} (173) 

with the usual definition of averages, such as the one given in Eq. ()154p . 

Generally for a field </>(x) with renormalized mass m and correlation length ^ = m~^, wave 
function renormalization constant Z, and (Euclidean) propagator 

< Hx)m > = / T^f^ e-^'-' > (174) 

J (ZTTj" + 

one has for $ = 4>{x) 

<$^>= / < <j){x)<j){y) > = V f <^{x)^{0) > = V = V Zf . (175) 

Jx,y Jx ^ 

® It should be noted that in the case of the lattice Wheeler-DeWitt equation of Eqs. I|20|l and H2ip and, generally, 
in any lattice Hamiltonian continuous-time formulation, the lattice continuum limit along the time direction has 
already been taken. This is due to the fact that one can view the resulting 2-1-1 theory as originating from one where 
there exist initially two lattices spacings, at and a. The first one is relevant for the time direction and the second one 
for the spatial directions. In the present lattice formulation the limit at ^ has already been taken; the only limit 
left is a — >■ 0, which requires the existence of an ultraviolet fixed point of the renormalization group. 
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Thus the field fluctuation probes the propagator at zero momentum, which in turn is directly 
related to the renormaUzed mass (and thus ^) for the field in question. If the field ^ acquires a 
nonzero expectation value, the above result is modified to 



(176) 



involving instead the connected propagator. In the gravity case the quantity Atot plays the role of 
if the fiuctuation diverges (^ — )• oo) then one has a phase transition or an ultraviolet fixed point 

in quantum field theory language [T71 [30] . 

Without the curvature term in the Wheeler-DeWitt equation [/3 = for the wave function ^ in 

Eq. (|16ip ]. the area fiuctuation does not diverge, even when n is large and is simply proportional 

to g^. In this case one finds 



xaW = 0) 
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2n-l 7r2(4n- l)r(4n-2) 
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.2 , /o^ A (177) 



47r n 



1 2i6"-i3(2n + i)r(n)8 

Note the spurious singularity for the special case of the tetrahedron, n = 1/2. When the curvature 
term is taken into account one finds, from the full wave function ^ in Eq. (jl6ip and in the limit of 
large n, 

/ 9\ n2 / 9 1 

(178) 
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Note that the fiuctuation now appears to diverge as 5 — )• (see also Fig. [7]). Furthermore, XA is 
nonanalytic in the original Newton's coupling G = g^, which suggests that perturbation theory in 
G is useless. A divergence of the fluctuations as — )• implies that in this limit the correlation 
length diverges in lattice units, signaling the emergence of a massless excitation. 

Just as for the case of the average curvature [Eq. (I169p ]. an exact sum rule can be derived in 
the (Euclidean) lattice path integral formulation, relating the local volume fluctuations to the local 
curvature fluctuations. In the three-dimensional Euclidean path integral theory the following exact 
identity holds for the fluctuations [30] 
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> - < 



h >^ 
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2iVi = , (179) 



> - < J^'^y^ > 

h h 

where A^i is the number of edges on the lattice (further exact sum rules can be derived by considering 
even higher derivatives of the free energy In Zl with respect to the parameters Aq and k). Since the 
last equation relates the fluctuation in the curvature to fluctuations in the volumes, it also implies 
a relationship between their singular (divergent) parts. 



^ We noted previously that in our Hamiltonian formulation, the lattice continuum limit along the time direction 
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Figure 7: Area fluctuation XA g = v G for the octahedron and icosahedron, computed from 
Eq. (jl72p . Note the divergence for small g. 



According to the sum rule of Eq. (jl79p a divergence in the curvature fluctuation 

XR -<{Y, hih? > - < E '^^^^ (180) 

h h 

for the three-dimensional (Euclidean) theory generally implies a corresponding divergence in the 
volume fluctuation 

~ < ( Vhf > - < >2 (181) 

h h 

for the same theory. In our case a divergence is expected in 2 + 1 dimensions of the form 



XA ~ \g - gd 



(182) 



with exponent a = l — 5 = 2 — Si^, where 5 is the universal curvature exponent defined previously 
in Eq. (jl70p . and z/ the correlation length exponent. The latter is defined in the usual way [31} [32] 
through 



C ~ \9- gc 



(183) 



where ^ is the invariant gravitational correlation length. The scaling relations among various 
exponents (i^, 5, a) are rather immediate consequences of the scaling assumption for the singular 



has already been taken. This results in two lattice spacings, one for the time and one for the space directions, denoted 
here respectively by at and a, with the first lattice spacing aheady sent to zero. It is then relatively straightforward 
to relate volumes between the two formulations, such as V ~ atA. Relating curvatures (for example, ^i? in the 2 + 1 
theory vs the Ricci scalar R in the original three-dimensional theory) in the two formulations is obviously less easy, 
due to the presence of derivatives along the time direction. 
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part of the free energy, Fging ~ the vicinity of a critical point (for more detailed discussion see, 

for example, [191 [311 [32] ) . The preceding argument then implies, via scaling, that a determination 
of a provides a direct estimate for the correlation length exponent v defined in Eq. ()183p . Note 
that based on the results so far one would be inclined to conclude that for 2 + 1 gravity the critical 
point — as n — )• oo. Equation (|182p can then be rewritten either as 



a/u 

l/d 



XA ~ e'" (184) 



or, in a finite volume with linear lattice dimensions L ~ Nq ~ V^A. ~ (since A^a = 4n + 2), 
as 

since, for a very large box and g very close to the critical point Qc, the correlation length saturates 
to its maximum value ^ ~ L. Hence the volume- or n-dependence of x provides a clear and direct 
way to estimate the critical correlation length exponent u defined in Eq. (jl83p . 



9 Results for Arbitrary Euler Characteristic x 



The results of the previous sections refer to regular triangulations of the sphere (x = 2) and the 
torus (x = 0) in 2 + 1 dimensions. It would seem that one has enough information at this point 
to reconstruct the same type of answers for arbitrary x- In particular one has for the parameter /3 
[see Eqs. (08]) and dnU] 



/3 



27rx 



(186) 



relevant for the wave functions in Eqs. (jl35p or (jl47p . For the average total area one then finds, 
using the wave function expansion in Eq. (|16ip . 



< Atot >= g < 



2i-2"r(n-i)r(2n + i) 



(n + l)7rx 



r(n)3 



+ 



■ 4l-^"r(n-i)r(n+l)r(2n+l)^ 
^ n2r(n)6 



g^ (2n + 1) 



In the large n limit one obtains for the average area of a single triangle 



< Aa> 



a 



2^/27 



n 



16n 



+ 



(7r-2)x 
g"^ n 



1 + 



1 



4n(7r - 2) 



+ ... 



and for the average total area 



/2n 4(^-2)x 
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+ 



(187) 



(188) 



(189) 
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For the area fluctuation defined in Eq. ()173p one finds in the same large n limit 



XA = (l--) ^ + o(-) + (4-7r)./^^ + ••• . (190) 
\ 7r/4 \^ / \ nn g 

Again note that the fiuctuation appears to diverge as — )• 0, which implies that this is the more 
interesting limit, so from now on we will focus specifically on this limit. It is clear from the analytic 
expression for < Atot > in Eqs. (jl87p or (jlSSp that as n — )• oo, the gravitational coupling g{n), to 
this order in the Bessel expansion, has to scale like 

g{n) ~ ^ , (191) 



so that the expression for < A^ot > scales like n or N^, with the expression for < tIa > staying 
finite. 

The result of Eq. (jl90p for XA then implies 

XA ^ ~ n° (192) 

in the same limit n — t- cxd. In view of Eqs. ()187p and ()185p with n ~ A^a ~ L'^, this would imply 
2/i> — 3 = 0, and thus for the universal critical exponent u itself z/ = | = 0.666 to first order 
(m = 1) in the Bessel function expansion of Eq. ()16ip and = = 0.588 to the next order 
(m = 2) in the same expansion. 

With some additional work one can, in fact, completely determine the asymptotic behavior of 
various averages for large f3 (small ^f) and large n. First one notes that when m Bessel functions 
are included in the expansion for the wave function given in Eq. ()16ip . beyond the leading order 
one at strong coupling, one obtains a wave function which contains powers of /3 up to /3™. For a 
given fixed m one then finds for the average area per triangle the following asymptotic result 

< > ~ ^ ^ ^ , (193) 

g.im—L j2 2 

up to terms which contain higher powers of 1/n (making these less relevant in the limit n — )■ oo), 
and also up to terms which are less singular in g for small g. The requirement that the average 
area per triangle be finite as n — )• oo then requires that the coupling g itself should scale with n 
according to 

g{n) ~ . (194) 

^2(3m-l) 

For the area fluctuation itself one then computes in the same limit 

XA ~ ^ , (195) 
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again to leading order in 1/n and The requirement that g{n) scale according to Eq. ()194p then 
implies from Eq. (jl95p that the area fluctuation diverges in the limit n — )• 00 as 

m— 1 

XA{n) ~ nsm-i . (196) 

By comparing with Eqs. (jl84p and (jl85p one obtains immediately for the exponent 

a 2m -2 



V 3m — 1 

and, therefore, from the scaling relation q = 2 — 3 finally 

6m — 2 



(197) 



(198) 



11 m — 5 

One can now take the limit m — t- 00 [infinite number of Bessel functions retained in the expansion of 
Eq. (jl6ip ]. which leads to the exact result for the correlation length exponent of 2 + 1 dimensional 
quantum gravity 

V = — = 0.5454... . (199) 

The derivation shows that the exponent v does not seem to depend on the Euler characteristic x 
and, therefore, on the boundary conditions ^ Furthermore one can compare the above value for 
V with the (numerically exact) Euclidean three-dimensional quantum gravity result obtained over 
twenty years ago in |24) . namely v ~ 0.59(2). It would, of course, be of great interest to repeat the 
above Euclidean lattice calculation in order to refine the estimate and improve on the statistical and 
systematic uncertainty. The exponent v is expected to represent a universal quantity, independent 
of short-distance regularization details and, therefore, characteristic of gravity's universal scaling 
properties on distances much larger than the lattice cutoff. As such, it should apply equally to both 
the Lorentzian and the Euclidean formulation, and our results are consistent with this conclusion. 
Moreover, in 3 -|- 1 dimensions the exponent is a key physical quantity as it determines the power 
for the running of the gravitational constant G [33] and for the Euclidean theory it is known [30J 
that the universal scaling exponent is consistent with v = 1/3. 

It is perhaps worthwhile at this point to compare with other attempts at determining the critical 
exponent v in three-dimensional gravity. The latest and best results for quantum gravity in the 
perturbative diagrammatic 2 + e continuum expansion using the background field method [35^ [36] 
give in d = 3 (e = 1 and central charge c = 1) 

z^~^ = 1 + 7 + . . . (200) 
5 



*One might wonder if the value for v is affected by the choice of normalization in Eqs. H56|) and (|155p . It is easy to 
check that at least the inclusion of a weight factor y4'", with m integer, does not change the result given in Eq. (|199|) . 
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to two-loop order and, therefore, v ~ 0.625, with a substantial uncertainty of about fifty percent 
(which can be estimated for example by comparing the one- and two-loop results). On the other 
hand, truncated renormalization group calculations for gravity directly in three dimensions |371 [38] 
give to lowest order in the truncation (i.e. with the inclusion of the cosmological and Einstein- 
Hilbert terms only) the estimate 

and, therefore, in d = 3 the value v ~ 0.833. This last result is also affected by a rather substantial 
uncertainty (again as much as fifty percent), which can be estimated, for example, by including 
curvature-squared terms in the truncated expansion. Nevertheless, and in light of the uncertainties 
associated with the various methods, it is very encouraging to note that widely different calculations 
(on the lattice and in the continuum) give values for the universal scaling exponent v that are 
roughly in the same ballpark. 

From Eq. (jl99p one obtains the fractal dimension for a gravitational path in 2 -|- 1 dimensions 

v-^ = dp = — = 1.8333... (202) 

This is slightly smaller than the value for a free scalar field dp = 2, corresponding to the Brownian 
motion (or Wiener path) value. It is closer to the value expected for a dilute branched polymer in 
the same dimension \39\ 140). and the best match at this point seems to be the 0{n) vector model 
for n = —1. The exact value v = 6/11 for 2-1-1 gravity would then suggest a connection between 
the ground state properties of quantum gravity and the geometry of dilute branched polymers in 
the same dimension. 

In light of the results obtained so far it is possible to make a number of additional observations. 
First, note from Eq. (I194p that as n — )• oo, the critical point (or renormalization group ultraviolet 
fixed point) moves to g = 

g{n) ~ (203) 

For comparison, a variational calculation based on correlated product (Slater-Jastrow) wave func- 
tions [T] in 2 -|- 1 dimensions gave 

a! = '"^ , (204) 

where do > 2 was a parameter associated there with the choice of functional measure over edges. 
The variational result of Eq. (j204p can be compared directly with the result of Eqs. ()165p and ()203p . 
for X = 2 and = 2n + 2. Thus in both treatments the limiting value for the critical point for g 
in 2 -|- 1 dimensions is zero, (^c — ^ as the number of triangles — )• oo. 
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Physically, this last result implies that there is no weak coupling phase {g < gc, or in terms 
of Newton's constant G < Gc)- the only surviving phase for gravity in three dimensions is the 
strongly coupled one {g > gc or G > Gc)- Furthermore, the correlation length ^ of Eq. (jl83p is 
finite for g > and diverges at 5 = 0. In particular, the weak field expansion, which assumes g 
small, is expected to have zero radius of convergence, In a sense this is a welcome result, as in 
the Euclidean theory the weak coupling phase was found to be pathological and thus physically 
unacceptable in both three [23] and four dimensions [171 130| . It would seem, therefore, that the 
Euclidean and Lorentzian lattice results are ultimately completely consistent: quantum gravity 
in 2 + 1 dimensions always resides in the strong coupling, gravitational antiscreening phase; the 
weak coupling, gravitational screening phase is physically excluded. In addition, the exact value 
for v determines, through standard renormalization group arguments, the scale dependence of the 
gravitational coupling in the vicinity of the ultraviolet fixed point |34) . 



10 Summary and Conclusions 

In this paper we have discussed the form of the gravitational wave function that arises as a solution 
of the lattice Wheeler-DeWitt equation [Eqs. (|20p . (j2ip and (|34p ] for finite lattices. The main result 
was the wave function ^ given in Eqs. (jl35p . (jl47p and (jl53p with strong coupling limit (curvature 
term absent) corresponding to the choice of parameter /3=0. 

To summarize, and for the purpose of the following discussion, the wave function ^ given in 
Eq. (|153p can be written in the most general form as 

* ~ e-*^iFi (a, b,2ix) (205) 

up to an overall normalization constant A^, and with parameters related to various geometric 

' These circumstances are perhaps unfamihar in the gravity context, but are nevertheless rather similar to what 
happens in gauge theories, including compact Quantum Electrodynamics in 2 + 1 dimensions |41) . There, the theory 
always resides in the strong coupling or disordered phase, with a finite correlation length which eventually diverges 
at zero charge. 

Specifically, the universal exponent ly is related to the behavior of the Callan-Symanzik beta function for Newton's 
constant G in the vicinity of the ultraviolet fixed point by /3'(G)|g=Gc = —l/i^- Integration of the renormalization 
group equations for G then determines the scale dependence of G'(pi) or G(n) in the vicinity of the ultraviolet fixed 
point. Concretely, determines the exponent in the running of G. One finds G(n) ~ (C^Q) with □ = g^^VfiV^ 

the covariant d'Alembertian and ^ the renormalization group invariant correlation length. A broader discussion of 
renormalization group methods as they apply to quantum gravity can be found for example in [19) . 
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invariants 



X = 



I 



a 



2V2AG 



b 



V2X 



Atot 



V2X 



I 



(206) 



X 



G 



G 



In the above definitions one can trade, if one so desires, the total number of triangles A'^a for the 
total area 



Use has been made of the relationship between various coupling constants {g, G, (3, A, A) to reexpress 
the wave function -0 in slightly more general terms, as a function of the original couplings A and G 
appearing in the original form of the Wheeler-DeWitt equation [see for example Eqs. (j42p and (j44p ]. 
We did show that an equivalent form for the wave function ^ can be given in terms of Coulomb 
wave functions [see Eq. (jl47p ]. with argument 



and X defined as in Eq. ()206p . 

The above wave function is exact in the limit of large areas and completely independent of the 
weak field expansion. Nevertheless it is only correct to some low order in the same expansion in the 
limit of small areas. This situation was interpreted as follows. For large areas one has a very large 
number of triangles, and the short distance details of the lattice setup play a vanishingly small role 
in this limit. One recognizes this limit as being relevant for universal scaling properties, including 
critical exponents. For small areas on the other hand a certain sensitivity to the short distance 
properties of the lattice regularization persists, and thus a universal behavior is, not unexpectedly, 
hard to achieve. In any case this last limit, in the absence of a truly fundamental and explicit 
microscopic theory, is always expected to be affected by short distance details of the regularization, 
no matter what its ultimate nature might be (a lattice of some sort, dimensional regularization, or 
an invariant continuum momentum cutoff, etc.) 

In principle any well-defined diffeomorphism-invariant average can be computed using the above 
wave functions. This will involve at some point the evaluation of a vacuum expectation value of 
some operator 0{g) 




(207) 




(208) 




(209) 
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where d^[g] is the appropriate functional measure over three- metrics gij. Evaluating such an average 
is, in general, non-trivial, as it requires the computation of a (Euclidean) lattice path integral in 
one dimension less 

md{gm = M j d^Ag] d{g,,) exp{-5e//[5]} (210) 

with Seff[g] = — In and J\f a normalization constant. The operator 0{g) itself can be 

local, or nonlocal as in the case of the gravitational Wilson loop discussed in [l2]. Note that the 
statistical weights have many zeros corresponding to the nodes of the wave function ^, and that 
Seff is infinite there. 

In the previous sections we have shown that the wave function allows one to calculate a number 
of useful and physically relevant averages and fluctuations, which were later extrapolated to the 
infinite volume limit of inflnitely many triangles. It was found that these diffeomorphism- invariant 
observables point in 2 + 1 dimensions to the existence of a fixed point (a phase transition in statis- 
tical field theory language) at the origin, Gc = 0. One concludes, therefore, that the weak coupling 
(or gravitational screening) phase has completely disappeared in the lattice nonperturbative formu- 
lation and that the theory resides in the strong coupling phase only. By contrast, in the Euclidean 
theory it was found in [23] that the weak coupling or gravitational screening phase exists but is 
pathological, corresponding to a degenerate branched polymer. A similar set of results is found in 
the four-dimensional Euclidean theory, where the weak coupling, gravitational screening phase also 
describes a branched polymer. 

The calculations presented in this paper and in [Ij can be regarded, therefore, as consistent 
with the conclusions reached earlier from the Euclidean framework, and no new surprises arise 
when considering the Lorentzian 2 + 1 theory. Furthermore, we have emphasized before that the 
results obtained point at a nonanalyticity in the coupling at G = 0, signaling a strong vacuum 
instability of quantum gravitation in this dimension. In view of these results it is therefore not 
surprising that calculations that rely on the weak field, semiclassical or small G expansion run into 
serious trouble and uncontrollable divergences very early on. Such an expansion does not seem to 
exist if the non-perturbative lattice results presented here are taken seriously. The correct physical 
vacuum apparently cannot in any way be obtained as a small perturbation of fiat or near-flat 
spacetime. 

Let us add here a few final comments aimed at placing the present work in the context of 

previous calculations for the same theory. A number of attempts have been made over the years 

The nature of solutions to the lattice Wheeler-DeWitt equation in 3 + 1 dimensions will be discussed in a separate 
publication [43] . 
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to obtain an estimate for the gravitational wave functional ^[g] in the absence of sources. These 
generally have relied on the weak field expansion in the continuum, as originally done in [HI [9]. 
Thus, for example, one finds in 3 + 1 dimensions 

^e.p{-j//.^*F«^2^-«}. (211) 

where /i^^(k) is the Fourier amplitude of transverse-traceless modes for the linearized gravitational 
field in four dimensions. The above wave functional describes a collection of harmonic oscillator 
wave functions, one for each of the infinitely many physical modes of the linearized gravitational 
field. As in the case of the electromagnetic field, the ground state wave functional can be expressed 
equivalently in terms of first derivatives of the field potentials (the corresponding B field for gravity), 
without having to mention Fourier amplitudes, as 



TTi 



Clearly both of the above expressions represent only the leading term in an expansion involving 
infinitely many terms in the metric fiuctuation hij (and since they apply to an expansion about flat 
space, the cosmological constant term does not appear either). Now, in 2 + 1 dimensions the above 
expressions become meaningless, since there cannot be any transverse-traceless modes. The only 
expectation that remains true is that the wave functional should still depend on physical degrees 
of freedom only: it should be a function of the intrinsic geometry of 3-space and should not change 
under a general coordinate change. 

If one restricts oneself to local terms a number of invariants are possible in 2 + 1 dimensions. 
In principle, the wave function could depend on, besides the total area 

Atot = jcfx^ (213) 

4vrx = Jcfx^R, (214) 

rn = J (fx^R^ (215) 

with n an integer. The latter result follows from the fact that in d = 2, both the Riemann and 
Ricci tensors only have one component, related to the scalar curvature, 

R/iupa — 2^^9fMcr9up 9pp 9iycr) i Rpv — 2 ^9p,u • (^-'-^) 



and curvature 



other invariants such as 
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Nonlocal terms are possible as well involving inverse powers of the covariant d'Alembertian □, but 
these do not seem to play a significant role in the lattice theory. 

Now, the relevant Euclidean theory for the present work is, of course, gravity in three (2 + 1) 
dimensions. But in three dimensions the Riemann and Ricci tensor have the same number of 
algebraically independent components (6) and are related to each other by 

R'\. = e^^'' eA.p [R". - \ (217) 

The field equations then imply, using Eq. (j217p , that the Riemann tensor is completely determined 
by the matter distribution implicit in T^y^ 

+ g^a Tup - gup T^a + T {Qp^ Qup - Qpp gua)] (218) 

In empty space Tpu = 0, which then implies for zero cosmological constant the vanishing of Riemann 
there 

Rpypa = . (219) 

As a result in three dimensions classical spacetime is locally flat everywhere outside a source, gravi- 
tational fields do not propagate outside matter, and two bodies cannot experience any gravitational 
force: they move uniformly on straight lines. There cannot be any gravitational waves either: the 
Weyl tensor, which carries information about gravitational fields not determined locally by matter, 
vanishes identically in three dimensions. 

What seems rather puzzling at first is that the Newtonian theory seems to make perfect sense 
in d = 3. This can only mean that the Newtonian theory is not recovered in the weak field limit of 
the relativistic theory. To see this explicitly, it is sufficient to consider the trace-reversed form of 
the field equations, 

Rpu = SttG (^Tpu - gpuT^ (220) 

with T = T\, in the weak field limit. In the linearized theory, with hpu = g^u — ^Ipu, and in the 
gauge V \h^p — ^Vph^x — 0> ™^ obtains the wave equation 

Uhpu = -levrG {jpu - (221) 

with Tpu the linearized stress tensor. After neglecting the spatial components of Tpp in comparison 
to the mass density roo, and assuming that the fields are quasistatic, one obtains a Poisson equation 
for /loo, 

V2/ioo = -levrG^^Too (222) 
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In four dimensions this is equivalent to Poisson's equation for the Newtonian theory when one 
identifies the metric with the Newtonian field in the usual way via /iqo = But in three 

dimensions such a correspondence is obstructed by the fact that, from Eq. (j222p . the nonrelativistic 
Newtonian coupling appearing in Poisson's equation is given by 

GNewton = ^ G (223) 

and the mass density roo completely decouples from the gravitational field /iqo- As a result, the 
linearized theory in three dimensions fails to reproduce the Newtonian theory. 

In a complementary way one can show that gravitational waves are not possible either in the 
linearized theory in three dimensions. Indeed the counting of physical degrees of freedom for the 
d-dimensional theory goes as follows. The metric g^y has + 1) independent components; the 
Bianchi identity and the coordinate conditions reduce this number to \d{d + 1) — d — d = 3), 
which gives indeed the correct number of physical degrees of freedom (two) corresponding to a 
massless spin two particle in d = 4, and no physical degrees of freedom in d = 3 (and minus one 
degree of freedom in d = 2, which is in fact incorrect). Nevertheless, investigations of quantum 
two-dimensional gravity have uncovered the fact that there can be surviving degrees of freedom in 
the quantum theory, at least in two dimensions. The usual treatment of two-dimensional gravity 
|44| starts from the metric in the conformal gauge g^u{x) = e'l'{x)g^y, where g^j_y is a reference 
metric, usually taken to be the flat one. The conformal gauge-fixing then implies a nontrivial 
Faddeev-Popov determinant, which, when exponentiated, results in an effective Liouville action, 
with a potential term coming from the cosmological constant contribution. One would, therefore, 
conclude from this example that gravity in the functional integral representation needs a careful 
treatment of the conformal degree of freedom, since in general its dynamics cannot be assumed 
to be trivial. The calculations presented in this paper show that this is indeed the case in 2 -1- 1 
dimensions as well. 
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